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The dynamics of a two-dimensional vortex are analyzed within the framework of the nonlinear 
Schrodinger equation. Both a bare vortex and a vortex with an external mass trapped in a finite- 
sized core are considered. The bare vortex motion is found to be damped at all frequencies, while 
the finite core has a single resonant frequency. The force exerted by the fluid on the finite core can 
be expressed as a sum of dissipative and Magnus forces for sufficiently low frequencies, even when 
the core is small. 
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I. INTRODUCTION 



Vortices represent an important class of excitations in 
many-body systems. They are associated with quantized 
circulation and dissipative flow in superfluids flux 
penetration and the breakdown of superconductivity in 
type-II superconductors |J, and the Kosterlitz-Thouless 
phase transition in two-dimensional systems 0. The 
study of vortices in classical fluids is also well-developed, 
and relevant material can be found in any good hydro- 
dynamics textbook. Finally, the interested reader may 
wish to consult Ref. [0], a recent monograph covering 
vortices, vortex lines, and vortex dynamics in all manner 
of nonlinear media. 

The aim of this paper is to address two questions per- 
taining to vortex dynamics in a two-dimensional boson 
fluid. First, we would like to determine whether there 
is an undamped or weakly damped mode with a moving 
vortex core; here the core velocity is measured relative to 
the fluid at infinity. Second, we want to understand the 
dynamics of the core when an external mass is trapped in 
it. A commonly-used phenomenological model expresses 
the force on a moving vortex core as the sum of a damping 
force acting parallel to the core velocity (again, relative 
to distant fluid) and a Magnus force acting perpendicular 
to it H . This type of model can describe classical vortex 
dynamics; it is also experimentally known to be valid in 
superfluid 4 He, at least when the core is a macroscopic 
object jjj. We would like to know whether such a model 
also describes the dynamics at shorter length scales, on 
the order of the size of a bare vortex core. We work 
within the nonlinear Schrodinger approximation, which 
is applicable in the limit of a dilute, weakly interacting 
Bose gas. A recent publication by Demircan, Ao, and Niu 
addressed similar questions within this framework Q ; we 
found that this work contained several mathematical er- 
rors invalidating its main results. 

The outline of the paper is as follows. We first present 
the basic vortex model and derive equations of motion 
for the core and the phonon modes to which it couples. 
We then show how these equations can be modified to 



include a massive object in the core. Next we discuss the 
qualitative features of the normal modes of the system 
and use numerical methods to find them. Finally we 
present our results and conclusions. 



II. VORTEX MODEL 

We start with the following Lagrangian ||, which de- 
scribes a two-dimensional system of bosons interacting 
via a delta- function pseudopotential: 



L= d 



~dt ~ 2m 



a) 

Here m is the boson mass, fi the chemical potential, 
and A the interaction strength. This can be put in di- 
mensionless form by a rescaling of variables. We let 
W -> (aVA) 1 / 2 *, r -> (h 2 /mfi)^ 2 r, and t -> (h/n)t to 
obtain 
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The prefactor is irrelevant for the present (classical) anal- 
ysis, since the classical equations of motion are invari- 
ant under a rescaling of L. However, note that the ac- 
tion J L dt acquires a prefactor of h /m\, which is much 
larger than h provided that A <§; h 2 /m; this defines the 
weakly interacting limit. In this limit the correspondence 
principle applies, at least naively, and quantum fluctu- 
ations around the classical behavior are expected to be 
small. The Euler-Lagrange equation for the complex field 
*(r, t) is 
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This is the well-known nonlinear Schrodinger equa- 
tion, derived for the imperfect Bose gas by Gross and 
Pitaevskii jlO], [H]. It has since been applied to vortex 
lines in superfluids by a number of authors (see, e.g., 
Refs. [[l2| [H| Q). It can be given a hydrodynamic 
interpretation by making the Madelung transformation 
$ = y/pe^, where p and </> are the fluid density and ve- 
locity potential (v = V</>) respectively. The equations of 
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motion for these variables are 
dp 



dt 



V ■ (pv), 



dt 



(pv) = V <T, 



(4a) 
(4b) 



where d/dt = d/dt + v • V is the convective derivative 
and the components of the stress tensor er are 



Except for the derivative terms in the stress tensor, these 
are the equations of motion for an ideal classical fluid 
with pressure p(p) = ^(p 2 — 1). 

Equation ([?]) has time-independent vortex solutions of 
the form f{r)e m9 for any integer n. We consider only 
ri = 1, in which case the solution is 

*o(r) = f(ry e , (6) 

where f(r) satisfies 

fir) + I. f^) + U- 2/(r) 2 - 1) f(r) = 0. (7) 

The asymptotic behaviors of the solution for large and 
small r can be expressed as power series: 
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' H as r -> 0, 

as r — > oo. 

(8) 

The value of A = f'(0) is not determined by the asymp- 
totics and must be found numerically. We used a shoot- 
ing method and found A ~ 0.825, which does not agree 
with the value of \[2 given in Rcf. ||, but does agree 
with the numerical work of Kawatra and Pathria [pL4j] . 

In order to investigate vortex motion, we consider per- 
turbations around a single- vortex solution with a moving 
core. The new dynamical variables are defined by 



tf(r,t)=*o(r-r (t)) + SV(r - r (i), t 



(9) 



where ro(t) is the location of the vortex core. Changing 
variables in Eq. (|^) and linearizing in 8^> and dro/dt gives 
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We express iJ't as a sum over cylindrical harmonics: 
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It is also convenient to rep- 



resent r by a complex number r + = r • (x + iy). In 
Eq. ([To]) , only the m = and m = 2 modes couple to the 
dynamics of Tq. The other modes can be set to zero, and 



the evolution of the relevant modes is given by 
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where the primes represent differentiation with respect 
to r. We impose the constraint 5^(0, t) = 0, so that 
ro marks the true location of the vortex core. These 
equations then determine the motion of both the fluid 
and the core. In particular, since t he S ^ m and their 
time derivatives vanish as r — > 0, Eq. ( 11a ) implies that 



dr+ 
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III. FINITE CORE 



We also would like to introduce an additional mass M, 
coupled to an arbitrary external force F(t) and confined 
to the center of the vortex. The question arises as to how 
this mass (which might physically represent a foreign par- 
ticle trapped in the core) should interact with the fluid. 
We considered using a point mass constrained to lie at 
the vortex core, but found this model to be internally 
inconsistent for the following reason. If the perturbation 
is bounded, as the validity of our linearized approach 
requires, the fluid exerts no force on a point mass. (This 
is justified in the next paragraph.) The equation of mo- 
tion for a point mass is therefore Mtq = F. On the other 
hand, for the mass to remain at the vortex core, Eq. Jl2] ) 
needs to be satisfied as well. The two equations for ro 
cannot be satisfied simultaneously for arbitrary F, so the 
model is ill-defined. Instead, we will treat the core as 
a hard disk of finite radius a. This is a simple approxi- 
mation to a more realistic interaction between the fluid 
and a trapped particle, such as a Lennard- Jones poten- 
tial. The fluid density must vanish inside the core, so 
the stationary vortex solutions now satisfy /(a) = 0, and 
(^(r, t) = for r — a. Figure |l| shows f(r) for the three 
cases we considered: the bare core, a small core (a = 0.2), 
and a large core (a = 2.0). The time dependence of small 
perturbations in the fluid is still given by Eqs. (p"l"|). 

For momentum to be locally conserved, the force on the 
core due to the fluid must equal the momentum crossing 
the boundary line r = a per unit time. This condition 
can be expressed in terms of the stress tensor: 
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FIG. 1: Square root of the fluid density for the bare core 
(solid curve), small core (dashed curve), and large core (dot- 
ted curve). 

where n is the boundary's outward unit normal and ds 
is the line element along the boundary. Integrating this 
expression around the circumference of the core and using 
the form of the stress tensor given in Eq. (|J) , the force 
exerted on the mass by the fluid is found to be 

F fluid = F fluid • (x + iy) 

= 7raf(a){s%(a) + (S^Y(a)). (14) 

Note that as a — > 0, this force vanishes unless (5*o + (5*2 
I 



the bare core will be considered as a special case with 
a = M = F + = 0. Equations (p~6|b-c) are coupled in- 
homogeneous linear ODEs for S^oyr) and (5*2(?") which 
must be solved subject to homogeneous boundary condi- 
tions: the (5* m must remain finite as r — > oo, and must 
vanish as r — > a. The most general solution will have 
the form of a particular solution to the inhomogeneous 
equations plus a linear combination of solutions to the 
corresponding homogeneous equations (obtained by set- 
ting r + = 0). 

We now consider the solutions to the homogeneous 
equations. Given the order of the equations, there must 
be four such solutions. As r — > their asymptotic behav- 



diverges at least logarithmically, as asserted above. Here 
we mention that Ref. || draws its conclusions from the 
point-mass model, which we have argued to be inconsis- 
tent. Their presentation of this model differs from ours 
in that they introduce a mass term directly into the La- 
grangian. This hides the inconsistency but does not re- 
move it: it leads to a seemingly nontrivial equation of 
motion for ro which is actually equivalent to Mtq = 0. 
A subsequent error (in joining inner and outer asymp- 
totic solutions) disguises this equivalence, enabling the 
derivation of quantitative results that are essentially un- 
supported. 

The equation of motion for the core mass is now 

Mr + = F« uid + F+, (15) 

where the external force has been written as a complex 
number F+ = F • (x + iy) . Since the equations of motion 
(Eqs. ([□]) and are linear in F+, r+, <5* , and (5*2, 

we will seek solutions in which all of these variables are 
proportional to e lwt . The general solution can be then 
expressed as a sum over normal modes in the usual way. 



IV. NORMAL MODES 

For a finite core, the mode equations for frequency ui 
are 



(16a) 



(16b) 



(16c) 



iors are as follows: (a) (5*o ~ r 6 and (5*2 ~ const, x r 2 ; 
(b) (5* ~ 1 and (5*2 ~ const, x r 4 ; (c) (5* ~ mr 
and (5*2 ~ const, x r 4 lnr; and (d) (5*o ~ r 2 and 
(5*2 ~ const, x 1/r 2 . In the opposite limit, as r — > oo, 
two solutions are oscillatory, representing incoming and 
outgoing waves, one solution decays exponentially, and 
one solution grows exponentially. (The asymptotic solu- 
tions were found by Fetter jL2|.) All other points (r ^ 
or oo) are regular points of Eqs. (p^b-c), so the (5* m and 
their derivatives can be specified freely. 

In the case of the bare core, only inner solution (a) 
conforms to the boundary condition at the origin. If this 
function is followed to infinity, it will decompose into a 



= F + 



<5* ' + h% + (2 - 4/ 2 - 2lo) 5*o - 2/ 2 (5*^ = -cur+ (f + t 
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linear combination of the four outer solutions, and in 
general (for arbitrary uf) the coefficient of the exponen- 
tially growing solution will not vanish. In other words, if 
the inner boundary condition is satisfied then the outer 
one will not be. Therefore, with a bare core there is no 
nonzero solution to the homogeneous equations, and the 
solution to the inhomogeneous equations will be unique 
and proportional to r+. 

In the case of a finite core, two linearly independent 
solutions satisfy the boundary condition at r = a, since 
6^' a (a) and (S^^Yi ) are independent free parameters. 
Both have some coefficient of overlap with the growing 
solution, so a linear combination can be formed for which 
the overlap is zero and both boundary conditions are sat- 
isfied. The general solution to the inhomogeneous equa- 
tions is in this case the sum of a term proportional to 
r + and a term proportional to another free parameter. 
Physically we expect to find a unique solution where the 
mass is driven solely by the external force and the fluid 
contains only outgoing waves. This solution should be 
proportional to r + . The free parameter associated with 
the homogeneous solution can be taken to be the ampli- 
tude of incoming waves. Once the solutions are found, 
the force exerted on the moving mass by the fluid is given 
by Eq. ©. 

At this point we resort to numerical techniques to de- 
termine the solutions for the two cases. 



V. NUMERICAL METHODS AND RESULTS 

No single numerical technique proved capable of solv- 
ing Eqs. (|l6|b-c) over the entire range a < r < oo. 
Instead, a different method was used for each of three 
regions: a < r < a + A, a + 4 < r < 2ir/ui, and 
2it /lu < r < oo . In each region we found a single solution 
to the inhomogeneous equations and all four linearly in- 
dependent solutions to the homogeneous equations. We 
then chose the appropriate coefficients for the homoge- 
neous terms so that S^o and satisfied the boundary 
conditions and joined smoothly at the boundaries be- 
tween regions. 

For the innermost region we used direct numerical in- 
tegration, starting at r = a, to determine the homoge- 
neous solutions. Because of the exponentially growing 
solution, the integration is unstable and loses precision 
with increasing r. We cut off the integration at r = a + 4; 
the exact placement was somewhat arbitrary, but it can- 
not be much further from the core. For the outermost 
region we used an asymptotic analysis. All the outer so- 
lutions to the homogeneous equations have asymptotic 
expansions of the form 

(17) 

where k takes on one of the four values ± y2± 2Vl + oj 2 . 
Using the full outer asymptotic series for f(r), the first 



few terms of which are shown in Eq. (|8|) , the coefficients 
a n {k) and b n (k) can be found for arbitrarily large n. 
However, we found a\/ao and bx/bo to be of order l/u> 
for small u>, indicating that the n = terms in the series 
dominate for r ^> 1/co. We therefore kept only the n = 
terms and cut off the outermost region at r ~ 2ir/uj. For 
the middle region we put Eqs. (|l6|b-c) on a grid, discretiz- 
ing the derivatives in the most straightforward way. For 
small values of to this middle region becomes quite large, 
so the average grid spacing needs to be reasonably large 
as well, to limit the number of grid points. On the other 
hand, the solutions still vary rapidly near the boundary 
at r = a + 4, so the spacing must be much smaller there. 
A nonuniform grid was used to address both these prob- 
lems. 

Finally, note that an exact solution to the inhomoge- 
neous equations, for all three regions, is 

6*0 = \r + (/' + tj , 5% = \r + (f tj . (18) 

This solution has a straightforward physical interpreta- 
tion: it describes a stationary vortex centered at the ori- 
gin. To see this, we can refer back to Eq. (^), the defin- 
ing equation for 5^ . A stationary vortex is described by 
*(r, t) = * (r), so 

<y*(r-r (t),t) = *o(r)-*o(r-r (t)) 

= r (t)-V* (r) (19) 

to first order in i"o , and the angular components of this 
function are indeed given by Eq. (|l8|). 

We found that our procedure gave good results for a 
wide range of frequencies. The results were extremely 
insensitive to adjustment of the region boundaries, in- 
creased precision in the numerical integration, and the 
inclusion of more terms in the asymptotic expansions. 
Changing the number and spacing of grid points in the 
middle region, however, did cause small variations in the 
final results. In particular, for 10~ 3 < \u>\ < 5 we esti- 
mate our errors to be on the order of a few percent. 

First we analyzed elastic scattering. When incoming 
waves of unit power are scattered by a free core (F + = 0) , 
the core responds with circular motion of a particular am- 
plitude. Figure || shows the amplitude of response as a 
function of frequency for the bare, small, and large cores. 
Note that for a finite core the response depends on the 
mass; we show the results for M = only, but using 
M > causes no qualitative changes. A divergent re- 
sponse is seen only as u> — > 0. The response functions for 
the small and large cores each have a single zero. These 
correspond to resonances in which the driven motion (dis- 
cussed below) is undamped and the core motion is de- 
coupled from the phonon field. The response function 
for the bare core, in contrast, is everywhere nonzero; we 
conclude that there is no undamped or weakly damped 
finite- frequency mode. 

For the small and large cores, we also analyzed driven 
motion. When a finite core is driven by the external force 
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FIG. 2: Amplitude of circular motion in response to incoming 
waves of unit power and frequency ui. Results are for the bare 
core (solid curve), small core (dashed curve), and large core 
(dotted curve). Ln|r+| is shown. 
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FIG. 3: Power dissipated by the small core (solid curve) and 
large core (dotted curve) in driven circular motion with unit 
speed and frequency lo. LnP is shown. 



with unit speed, power is radiated as outgoing waves. 
The power is equal to that put into the system by the 
component of F parallel to the velocity: for r + = 1/\lo\, 
corresponding to circular motion with unit speed, it is 
given by P — F • dr /dt — ± lmF + , where the sign is 
+/— for positive/negative frequencies. This quantity is 
shown in Fig. I as a function of frequency. It is in- 
dependent of M, and so provides a better measure of 
damping than the elastic scattering response. For the 
small core, the radiated power vanishes at lo ~ 1.5. At 
this frequency there are neither incoming nor outgoing 
waves; the perturbation is localized and dies off expo- 
nentially with increasing r. There is a similar resonance 
for the large core at w ~ 0.13. For both cores, the mo- 
tion is undamped as lo — > 0. Results for other core sizes 
(not shown) suggest that these features are general: there 
is always a single positive-frequency resonance, the fre- 
quency of which varies inversely with a, and the damping 
always approaches zero for low frequencies. 

Figure |^ shows the component of the force exerted by 
the fluid on the mass in the —dro/dt x z direction. If the 
only contributor were the Magnus force, 



^Magnus 



= -2tt 



dr Q 
dt 



x z, 



(20) 



this component would be identically 27r, which is marked 
on the vertical axis. Because we are driving the core with 
unit speed, the dependence on lo should vanish. Instead 
we see that the Magnus effect dominates only in the limit 
as lo — > 0; for \lo\ > 1/10, it can no longer account for the 
perpendicular component of the force. 



VI. CONCLUSIONS 

From our data we can answer both of our original ques- 
tions. First, there is no undamped mode associated with 
the motion of a bare vortex. Second, the Magnus effect 
plays an important role in the dynamics of a finite core 




CO 

FIG. 4: Force in the —dro/dt x z direction exerted by the fluid 
on the small core (solid curve) and large core (dotted curve) 
in driven circular motion with unit speed and frequency u. 
The mark on the vertical axis corresponds to Fx = 2n. 



for the frequency range \lo\ < 1/10, even when the core is 
small. The Magnus force has been shown to be a general 
consequence of vortex motion in the adiabatic phase ap- 
proximation Jl5| . We point out that this approximation 
is not valid when applied to the nonlinear Schrodinger 
equation, at least for an infinite system. To see why it 
must fail, again consider the solutions to Eqs. (p^b-c) in 
the limit as r — > oo. The particular solution given by 



Eq. (18), as already stated, corresponds to a stationary 
vortex. The phase field of distant fluid does not follow the 
instantaneous location of the vortex, as it does within the 
adiabatic phase approximation, but rather maintains a 
fixed center. The most general solution to the mode equa- 
tions has the same property, as the solutions to the cor- 
responding homogeneous equations are either oscillatory 



or exponentially decaying, as discussed in Section IV, 
and cannot cancel this power-law behavior. The failure 
of the adiabatic phase approximation is also expected on 
physical grounds: since the low-frequency phonons travel 
with finite speed (the speed of sound c = 1 in our units), 
the fluid cannot respond to low-frequency disturbances 



6 



at a distance d in a time shorter than l/d. For this re- 
sponse time to be shorter than the time scale associated 
with the disturbance, we must have d -C 1/w; at larger 
distances the phase of the fluid can no longer follow the 
core. Nevertheless we recover the Magnus effect for low 
frequencies, in agreement with results dependent on the 
adiabatic phase approximation, suggesting that the dis- 
crepancy at large distances is not crucial. For high fre- 
quencies, the Magnus force ceases to play a role. Note 
that this transition occurs while the phonon wavelength 
is still very large compared to the size of the core: for 



lu = 1/10, the phonon wavelength is around 60. For 
superfluid 4 He, the transition frequency is 20-30 GHz, 
corresponding to a phonon wavelength of 100 A. 



Acknowledgments 

This work was performed with the kind support of Cor- 
nell University. 



[1] For the published version of this article, see: M. J. Quist, 

Phys. Rev. B 60, 4240 (1999). Copyright 1999 by the 

American Physical Society. 
[2] R. P. Feynman, in Progress in Low Temperature Physics, 

edited by C. J. Gorter (North-Holland, Amsterdam, 

1955), Vol. 1. 

[3] M. Tinkham, Introduction to Superconductivity (McGraw 

Hill, New York, 1996). 
[4] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 

(1973). 

[5] L. M. Pismen, Vortices in Nonlinear Fields (Clarendon, 
Oxford, 1999). 

[6] R. J. Donnelly, Quantized Vortices in Helium II (Cam- 
bridge University Press, Cambridge, 1991), Chap. 3. 
[7] W. F. Vinen, Proc. R. Soc. A 260, 218 (1961). 



[8] E. Demircan, P. Ao, and Q. Niu, Phys. Rev. B 54, 10 027 
(1996). 

[9] J. W. Negele and H. Orland, Quantum Many-Particle 
Systems (Addison- Wesley, Reading, 1987). 

[10] E. P. Gross, Nuovo Cimento 20, 454 (1961). 

[11] L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 640 (1961) [Sov. 
Phys. JETP 13, 451 (1961)]. 

[12] A. L. Fetter, Phys. Rev. 138, A709 (1965). 

[13] A. L. Fetter, Phys. Rev. 138, A429 (1965). 

[14] M. P. Kawatra and R. K. Pathria, Phys. Rev. 151, 132 
(1966). 

[15] P. Ao and D. J. Thouless, Phys. Rev. Lett. 70, 2158 
(1993). 



